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ABSTRACT 

Spatial variability of canopy temperature has been shown to provide an useftil criteria for irriga- 
tion scheduling, and a knowledge of this variability is needed for comparing infrared radiometric 
observations of differing spatial resolutions. A soil-plant-atmosphere model for com (Zea mays L.) 
together with the scaling theory for soil hydraulic heterogeniety are used to study the sensitivity 
of spatial variation of canopy temperature to field-averaged soil texture and crop rooting charac- 
teristics. The soil-plant-atmosphere model explicitly solves a continuity equation for water flux 
resulting from root water uptake, changes in plant water storage and transpirational flux. Dy- 
namical equations for root-zone soil water potential and the plant water storage models the pro- 
gressive drying of soil, and day-time dehydration and night-time hydration of the crop. The 
statistic of scaling parameter which describes the spatial variation of soil hydraulic conductivity 
and matric potential is assumed to be independent of soil texture class. The field-averaged soil 
hydraulic characteristic are chosen to be representative of loamy sand and clay loam soils. Two 
rooting characteristics are chosen, one shallow and the other deep rooted. The simulation shows 
that the range of canopy temperatures in the clayey soil is less than IK, but for the sandy soil the 
range is about 2.5 and 5. OK, respectively, for the shallow and deep rooted crops. Irrigation 
scheduling bas'd on a fixed canopy temperature variability would lead to variable depletion of 
available water, depending upon soil texture and crop rooting characteristics. Because of a smaller 
variability in clayey soils, one may use the ground-based radiometric observations to ‘calibrate’ 
coarse resolution observations (as would be available from a air-borne sensor), provided that a soil 
hydraulic heterogeniety is the only important heterogeniety within the resolution area. 
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SPATIAL VARIATION OF CORN CANOPY TEMPERATURE 
AS DEPENDENT UPON SOIL TEXTURE AND 
CROP ROOTING CHARACTERISTICS 


INTRODUCTION 

Production of agricultural crops in arid and semi-arid regions of the world is highly dependent 
upon irrigation. In humid regions irrigation is also being used to increase yield. While water used 
for irrigation is increasing, the water available for irrigation is decreasing (aquifer recharge being less 
than depletion). There is an increasing emphasis to maximize the benefit from each unit of irriga- 
tion water. The key question posed in irrigation scheduling is when and how much water is to be 
applied to offset the yield-reducing stress conditions of a crop. 

Recently Jackson (1982) has reviewed the irrigation scheduling techniques based on soil, 
meteorologic, and plant factors. He argues that while soil and meteoro logic factors used in irrigation 
scheduling stride to indicate crop stress and needed irrigation water, direct observation of some 
plant parameter should be a superior approach. The plant parameters which have been studied as 
indicators of plant stress are leaf water potential, stomatal resistance, and leaf or canopy tempera- 
tures. The leaf based measurements (potential, resistance and temperature) are time consuming 
when fields of the size of hectare or larger are to be sampled to obtain a representative value. 
Therefore, according to Jackson (1982) the most promising approach to irrigation scheduling would 
be the measurements of canopy temperature using infrared radiometers. A further highly appealing 
nature of irrigation scheduling or identifying crop stress using infrared radiometers is that the 
technique could be applied to large areas using airborne and spacebome infrared sensors. 

At high soil water potentials the canopy temperatures are independent of soil texture and the 
actual soil water potential; the canopy temperatures being determined by atmospheric evaporative 
demand and plant characteristics. Since with decreasing soil water potential the canopy tempera- 
ture increases. Tanner ( 1 963) argued that it may be possible to infer soil water status from the 
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canopy temperatures. In order to deviate the effects of diurnal and day-to-day variations in 
weather conditions on the canopy temperature, Aston and van Bevel (1972) addressed the need of 
a reference canopy temperature. The reference temperature could be of a nearby well-watered plot. 
Alternately, they suggested that a cropped field can be its own internal reference, because, all 
fields being heterogeneous, a differential soil drying would introduce a spatial crop temperature 
variability. The magnitude of this variability may signal the onset of crop water deficits. 

The problem of accounting the effect of variable weather conditions on a crop canopy temper- 
ature has been addressed further by Idso et at. (1981a) and Idso (1982), among others. UaottaL 
(1981a) observed that during a significant portion of dear sky daylight period the canopy temper 
a turn of well watered crops are determined largely by air and dew point temperatures. Idso (1982) 
gives empirically derived linear regression equations for calculating unstress canopy temperature of 
26 crops in terms of air and dew point temperatures. These equations could be used to infer the 
crop water status from infrared radiometric observations knowing concurrent air and dew point 
temperatures. (Choudhury, 1982 and 1983). In remote sensing applications, however, it may be 
difficult to acquire the ancillary weather data of concurrent air and dew point temperatures. 

The alternate approach to irrigation scheduling suggested by Aston and van Bavel ( 1 972) 
which would not require any weather data, namely the use of spatial variability of canopy tempera- 
tures, was recently field tested by Clawson and Blad (1982). In sprinkler irrigated com plots of 
Valentine fine sand underlying day horizons of dimensions 18.3 X 36.6 m, Clawson and Blad 
( 1 982) found that irrigation with CTV * 0.8 K (where CTV * maximum minus minimum observed 
canopy temperatures) reduced the amount of irrigation water by about 54% with only 5% reduction 
in yield as compared to a well-watered plot. This saving of irrigation water is certainly worth 
noting, and while a plot-averaged canopy temperature would indicate the average soil water status 
of the field, the field variability of the canopy temperature indicates the differential root water 
uptake. The distribution of root-zone soil water within a field would lead of a differential stress 


conditions ami hence a spatial variation of crop yield (Bresier et at, 1981 ). Irrigation using canopy 
temperature variability should be able to minimize the spatial variation of crop stress, and it follows 
from MitscherUch’s equation for plant growth (Ware etal, 1 982) that the total yield of a field 
could be increased. 

Apart from irrigation scheduling application, a knowledge of spatial variability of temperature 
is also needed in comparing radiometer observations of differing spatial resolutions (for example, 
ground-based and satellite-borne sensors), as has been pointed out by Vauclin etal. (1982). 

The degree of heterogeniety in agricultural fields is highly variable. The heterogeniety could be 
both deterministic (topographic) and stochastic, and at present there are no general framework for 
describing the variability of soil hydraulic properties. However, by sacrificing some generality, the 
scaling theory of Philip (1967) has provided a simple physically based model to describe a stochastic 
heterogeniety (Warrick etal.. 1977; Russo and Bresier, 1980). In this theory, the spatial variability 
of unsaturated hydraulic conductivity and matric potential is described by a frequency distribution 
and autocorrelation functions of a single stochastic variable, called the scaling parameter. The 
frequency distribution of the scaling parameter has been found to be normal (Russo and Bresier, 
1980) and log-normal (Warrick et al.. 1977), and the auto-correlation length varied from a few tens 
of meters (Vieira et al., 1981) to less than a meter (Sisson and Wierenga, 1981). Spatial variability 
of infiltration rates in a Yolo loam agricultural field studied by Vieira etal. (1 981 ) is quite signif- 
icant in showing that there were no correlation of the rates when the soil samples were separated 
by SO m or more, although the field was considered to be " exceptionally uniform " for agricultural 
purposes. Peck et al. (1977), Sharma and Luxmoore ( 1 979) and Lascano and van Bavel (1982), 
among others, have discussed the application of scaling theory to soil water flow problems. 

Assuming a scale heterogeneous medium with uncorrelated log-normal distribution of scaling 
parameter. Choudhury and Federer (1983) studied the effect of soil heterogeniety on a com canopy 
temperature. The simulation model of Choudhury and Federer solves the continuity equation of 
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water flux in the soil-plant-stmosphere system by relating the canopy reaistance and the rate of 

i 

change of plant tisaue water with the leaf water potential. The dynamical equations of soil- and 
plant tissue-water were integrated to simulate the effects of progre s s iv e soil drying and plant dehy- 
dration on the canopy temperature. For a field averaged soil hydraulic properties repreeentive of a 
sandy loam soil and a fixed rooting density, the simulation showed that during the first several days 
of soil drying the field averaged canopy temperature increases fairly slowly but the canopy tempera- 
ture variability shows a distinct onset with an almost exponential rise to a maximum value, after 
which it decreases. These simulated trends were found to be in qualitative agreement with the 
observation of Gardner et aL (1981) and Clawson and Blad (1982). 

The present simulation is of the nature of a sensitivity study of spatial variation of a com 
canopy temperature. The factors considered in this study are field-averaged soil texture and crop 
rooting characteristics. The spatial heterogeniety of the fields, as described by the statistic of 
scaling parameter, is assumed to be identical for all soil textures. The present simulation addresses 
the question - For fields which are statistically identical with regard to spatial heterogeniety but 
differ in mean hydraulic and crop rooting characteristics what would be the range of canopy 
temperatures? 

Clearly, the range of canopy temperatures would be directly re la vent in comparing radiometric 
observations of differing spatial resolutions. Differences in the canopy temperature which may 
result from radiometer look angle (cf., Jackson, 1982) is not addressed in this paper. 

From field experiments Clawson and Blad (1982) had found an efficient irrigation scheduling 
scheme for com with a canopy temperature variability (CTV) of 0.8 K, but they questioned whether 
this value would be applicable to other field conditions. The question of prescribing an efficient 
irrigation scheduling scheme based on the results of present simulation is not addressed because 
there are no unique mathematical criteria currently available to provide such a prescription. 
However, the field-averaged values of the depletion level for available water, the stomatal resistance 
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and the leaf -water potential, all of which have been used for irrigation scheduling (see Jackson, 
1982), will be discussed. 

SIMULATION METHODS AND PARAMETER SELECTION 
Soil-Plan t-A tmosphere Model 

The dynamics of soil-plant-atmosphere system is simulated using Choudhury and Federer’s 
(1983) model. The model is schematically illustrated in Figure 1, and a brief description follows. 


The continuity equation of water flux in the soil-plant-atmosphere system, namely 


AR n + c p />,(«;-e,)/r, j ♦,-» + »„-» 
L, [ A + -r (r^ + % )/r a ] R, + Rp Rq 


(where R n is net radiation, e* and e a are, respectively, the saturated vapor pressures at air and dew 
point temperatures, A is the slope of saturated vapor pressure with respect to the air temperature, 
7 is psychrometric constant, Ly is the latent heat vaporization, r a and r a are, respectively, the 
boundary-layer resistances for vapor and heat transfers and the other symbols axe defined in 
Figure 1) is made to be a self-consistent equation for determining the plant water potential (<//) by 
expressing the canopy resistance (r c ) in terms of ^ via the following empirical equation 


1.4 X 10 4 17S 

Fc “lai (*-* c ) 1 + S+ 10 


( 2 ) 


where LAI is the canopy leaf area index, S is the global insolation and is the critical plant water 
potential for stomatal closure. The eqn. (2) is derived from the observations of Uchijima (1976) 
and Reicosky and Lambert ( 1 978), and following the suggestion of Kramer ( 1 969) that the leaf 
water potential and global insolation are the most important factors determining the stomatal 
resistance. 
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The day-time plant dehydration and night-time recharge is modeled through an empirical 
relationship between the water stored in the plant (Q) (which is the volume of plant water above 
permanent wilting per unit surface area) and the corresponding potential ** 

* q --450 X f- (131 X + 84)(l-f) (3) 


where 



and Q 0 is the amount of water stored at maximum potential (i.e., ■ 0). 

The soil resistance (R,) calculated from root depth and soil hydraulic conductivity according 
to Feddes and Rijtema (1972), and the plant resists, ce (Rp) is calculated from root length per 
unit area according to Cowan (1965) taking into account the plant stress induced nonlinearities 
(Boyer, 1971; Nulsen and Thurtell, 1 978). Monteith’s (1981) equation is used to calculate the 
boundary-layer resistance for heat transfer (r a ) from wind speed and air temperature. 


The dynamics of soil water potential (*// s ) and the water stored in the plant (Q) are governed 
by the equations 

z eff C <*s> ■ “ “ 4+^ ’ K ( * s> (4) 

<» 

where Z e ff is root depth, K (^ s ) is the soil hydraulic conductivity and C (^ s ) is the soil capacitance. 
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The above soil-plant-atmosphere model applies to a vegetation-covered soil for which soil 
evaporation is neflifible. The input hourly meteorolofical parameters are global insolation, wind 
speed and air and dew point temperatures; the plant parameters are the root depth, the leaf area 
index and the root length per unit area; the soil parameters are hydraulic conductivity and capaci- 
tance, which will be elaborated in the next section. 

With the above input parameters, Equation ( 1 ) is solved for the plant water potential ^ by 
initializing \ji % and Q. After ^ is calculated, the dynamics of and Q follow from Equations (4) 
and (S). The model thus generates successive steady states of plant water potential. The left hand 
side of eqn. (1 ) gives the transpirations! flux, and the canopy temperature follows from the energy 
balance equation. Further details of the model together with model validation by comparison with 
observed plant water potentials and stomatal resistances can be found in Choudhury and Federer 
(1983). A representative set of validation results are shown in Figure 2. The weather data of 
Figure (2b) taken from Reicosky and Lambert (1978) are used in the present simulation. 


Soil Heterogeniety 

In the present simulation, the hydraulic conductivity K (^ s ) and soil capacitance C are 
expressed in terms of soil water potential using the parametric equations of Clapp and Horn- 
berger (1978) as: 


K (^ s ) ■ K Mt 


*s 


C(* s ) 


^sat f ^sat i 

*sat b *s 


(6) 

(7) 


where b is an empirical constant and 0 Mt , and are the saturation values of, respectively, 
volumetric moisture, potential and conductivity. 
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According to Dagan and Bretler (1979) although all parameters of Equations (6) and (7) 
(namely, b, 9 uV and K ut ) are stochastic for a heterogeneous soil, the randomness could be 
confined to only K at and ^ Mt using the scaling theory as: 


Kut'a'Kat 

(8) 


(9) 


where the scaling parameter a is a stochastic variable and l£ ttt and *», are the field-averaged 
values. 

The field-averaged soil hydraulic parameters chosen for the present simulation correspond to 
loamy-sand and clay-loam soils (Gapp and Homberger, 1978); for loamy sand: b ■ 4.38, 0 ttt * 

0.41 , * -0.1 m and K Jgt ■ 1 .5 X 1(T* m/s, and for clay loam: b * 8.5, 6^ ■ 0.476, * 

-0.63 m and R sat » 2 X 10"* m/s. 

Based on the observations of Warrick et al. (1977) for soil series Panoche and Pima it is 
assumed that with 95% probability the spatial variation of the conductivity and the potential can 
be described by taking the scaling parameter a between 0.2 and 2.5 , with o*l for the mean soil 
hydraulic characteristics. The simulation results are, therefore, obtained with a“0.2, 1.0 and 2.5. 
(Clearly, if the spatial heterogeniety of the fields are not identical then for the same 95% proba- 
bility the bounds of a, namely 0.2 and 2.5, would have to be modified.) This approach to studying 
the effects of spatial variability on soil water problems is along the line of Bresler et al. (1979) and 
Lascano and van Bavel (1982). A slightly different approach was followed by Choudhury and 
Federer (1983) who did the dynamic simulation by stochastically generating a sequence of a 
values by prescribing a frequency distribution (the crude Monte-Carlo approach). The present 
approach is computationally much less time consuming. 

Com Rooting Characteristics 

The rooting characteristics (depth and density) cf all crops are highly variable, «u.d the 
significant factors affecting these characteristics are growth stage, soil temperature, texture and 


strength, irrigation management and planting density (see, for example, Kramer, 1969; Milthorpe 
and Moorby, 1974). 

At the silking time for com grown u< sandy soils Robertson et cl. (1980) observed the root 
length per unit area (L A ) of 5.7 X 10* nf 1 with no irrigation and 1.9 X 10 4 m“ 1 with light in- 
frequent irrigation. In a silt loam soil (which had a consolidated layer below 0.7S m) Mengel and 
Barber (1974) observed L A values of 9 X 10* and 1 .25 X 10 4 nf 1 , respectively, after 60 and 95 
days of planting. Barber (1971) observed significant variation in com root morphology with 
cultivation practice. From the observed variation of rooting density at different soil depths, Grimes 
et al. (1975) concluded that 75% of these variations could be accounted for in terms of soil 
strength. Reicosky and Lambert (1978) found that 85 to 90% of the roots of a mature com crop 
were in the 0.31 m thick Varina sandy loam Al horizon which was overlying a consolidated layer. 
However, the data of Grimes et al. (1975) show that for r-uure crops in permeable (low strength) 
soils 85 to 90% of the roots are in the surface meter of soil. For isolated com plants grown in 
loamy fine sand Taylor and Klepper (1973) observed root lengths greater than 16 km. If one 
neglects the root-to-root growth limiting interactions operative under field conditions then a 
root length of 16 km would give an L A value of about 1 X 10 s m~‘ for a planting density of 
7 plants/m : . 

The intention of the above discussion was to provide some quantitative results regarding the 
observed variation in com rooting characteristics. To simulate the effect of differing rooting 
characteristics on the canopy temperature the combination of effective rooting depth Z e ff (m) and 
root length per unit area L A (m _1 ) (i.e., (Z e ff, L A )) are chosen to be (0.3, 5 X 10* ) and (1.2, 

1 .5 X 1 0 4 ). The chosen range is fairly large, but one should recognize that a larger variation could 
occur. 
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SIMULATION RESULTS 

According to Kramer (1969) the single most important plant stress indicator is the leaf-water 
potential. The importance of the potential as an indicator of stress steins from the fact that cell 
division and growth, which are central to the ultimate crop yield, are related directly to the poten- 
tial (Hsio, 1973). Clark and Hiler (1973) tested several plant measurements (the potential, stomatal 
resistance and infrared temperature) in order to select the most effective crop stress indicator; the 
selected indicator was the leaf-water potential. Ehling et al. (1968) and Brady et al. (1974), among 
others, observed that leaf- and soil-water potentials am related linearly to each other. Idso et al. 
(1981b) found the leaf water potential to depend upon the atmospheric vapor pressure deficit. 

The simulated field averaged mid-day leaf and soil water potentials are shown in Figure 3. The 
relationship between the two potentials is faiily linear, although the linearity is more for the clay 
loam sod as compared to the loamy sand soil. Both both soils a marked deviation f:om linearity is 
seen to occur at low soil water pc tentials In Equation ( I ) the soil water potential appears linearly 
in the numerator and non-linearly in the denominator through the dependence of R s on the hy- 
draulic conductivity. At high soil water potentials the plant resistance (R p ) is several orders of 
magnitude larger than the soil resistance (R $ ), and, therefore, the dependence of soil water poten- 
tial in Equation ( 1 ) is essentially linear at high potentials. As the soil dries, the soil resistance 
increases and begins to compete with the plant resistance. Thus, with decreasing soil water poten- 
tial, the Equation ( 1 ) becomes a non-linear function of the potential. 

Since the plant resistance is the controlling factor for root water uptake when the soil water 
potential is high, differences in crop rooting characteristics wjuld be expected to affect the leaf 
water potential, as is seen in Fig. 3. Under a high evaporative demand, a crop with a higher rooting 
density would be able to transpire more efficiently and, therefore, would have a higher leaf water 
potential. A factor of two difference in the potentials of a deeply rooted and a shallow rooted 
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crop is seen in Figure 3. Thus, although leaf and soil water potentials are related lineraly, one 
will not be able to infer the soil water potential directly by measuring the mid-day leaf water 
potential. 

In order to account for the effect of atmospheric variability on the leaf water potential, 
Choudhury and Federer (1983) had defined an index, called the plant water stress index, as 

PWSI * y~~ 0 (10) 

where is the leaf water potential in the absence of any soil water deficit (the value assumed is 
for the day 1 of simulation) and is the critical potential for stomatal closure. The PWSIs cal* 
culated from \j/ values in Fig. 3 are shown in Fig. 4, which indicates that this index is also relatively 
insensitive to the crop rooting densities. If the long-held wisdom of irrigation scheduling based on 
60 to 70% depletion of available water is to be maintained then PWSI should not be allowed to 
exceed a value of about 0.35. Thus, instead of using the \p values directly to infer the soil water 
potential, one can use the index PWSI to infer the potential with less ambiguity. 

The relationship between pre-sunrise leaf and soil water potentials are shown in Figure 5. 
While a 1 : 1 relationship exists for a deep root crop, the leaf water potential of a shallow rooted 
crop is seen to be generally lower than the soil water potential. This inequality in the potentials 
results from day-time plant dehydration which get gradually restored during the night. Under a 
high evaporative demand a shallow rooted crop gets dehydrated more because their roots are less 
capable of extracting soil water so as to meet the evaporative demand. Since the restoration of 
lost plant water continues throughout the night, the leaf water potential reaches the highest value 
just prior to sunrise. Boyer ( 1 968) had noted that the leaf water potenti il before sunrise could 
be 10 to 20 m lower than the soil water potential (see also Brady et al., 1974). The leaf water 
potential during the n r:.. need: to be several meters lower than the soil water potential if water 
is to be ’pulled’ to the leaf for ceil growth and division and other metabolic use. 


11 


I 


The simulated field -averaged stomatal resistances and soil water potentials at 1 300 hour are 
shown in Figure 6. At igh soil water potentials the stomatal resistances depend upon rooting 
characteristics and atmospheric evaporative demand. As the evaporative demand increases or the 
rooting density decrease* the stomatal resistance increases (cf.. Fig. 6). During the drying period 
of an irrigation cycle, the stomatal resistance increases initially at a slow rate but then fairly 
rapidly. In a shallow rooted crop, the rapid increase in the stomatal resistance occurs earlier than 
for a deep rooted crop because of the relatively quicker depletion of root-zone soil water in a 
shallow rooted crop. The dependence of stomatal resistance on soil water potential, as could be 
derived from this figure, agrees well with the observations of Brady et aL (1975). 

Since in situ observations of the root-zone soil water potential are subject to errors in measur- 
ing the root distribution, a measurement of stomatal resistance has been suggested to be an indirect 
way of inferring this potential (cf., Brady et al., 1975). The present simulation although shows that 
the stomatal resistances are affected by the soil water potential, it questions any unambiguous 
determination of the potential by measuring the resistance only. Figure 6 clearly shows such an 
ambiguity at high soil water potentials, since both resistances 1 20 s/m and 340 s/m correspond to a 
soil water potential of about -3 m. A measurement of pre-sunrise leaf water potential appears to 
be the most direct approach to measuring the soil water potential (c.f., Figure 5.) For irrigated corn 
crops Sumayao et al. ( 1 980) observed the stomatal resistances to be lower than 500 s/m. The pres- 
ent imulation would suggest that this value of the stomatal resistance for irrigated com may not be 
applicable to other locations having a different evaporative demand or having a different rooting 
characteristics. Thus, in irrigation scheduling using stomatal resistance measurements, one would 
need additional information about an average evaporative demand and the crop rooting characteristics. 

The simulated field-averaged canopy-air temperature differences (8T) and soil water potentials 
at 1 300 hour are shown in Fig. 7. Similar to the above results for the leaf water potentials and sto- 
matai resistances, the canopy temperatures are seen to be affected by the crop rooting density, in 
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addition to the soil water potential. The relationship of 6T with the soil water potential is approx* 
imately linear. The canopy temperatures of a deeply rooted crop are lower compared to a shallow 
rooted crop. 

The above simulation results pertained to field-averaged leaf water potentials and stomata! 
resistances. While both of these plant based measurements have been used to characterize soil 
water conditions of agricultural fields, they are more time consuming than the canopy temperature 
measurements using infrared radiometers. Baaed on the suggestion of Aston and van Bavel ( 1 972), 
Clawson and Blad (1982) found that the spatial variability of canopy temperatures would be a good 
indicator of crop water deficiencies. 

The simulated range of 1300 hour canopy temperatures (i.e.. maximum minus minimum tem- 
perature) are shown in Figures 8a and 8b, respectively, for clay loam and loamy sand soils. Each 
figure contains results for two rooting characteristics, and the depletion of 70 percent of available 
soil water is annotated. The depletion of 70 percent of available water has been used as a criterion 
for irrigation scheduling. 

The canopy temperature variability is clearly a function of soil texture and rooting charac- 
teristics. In coarse textured sandy soils the effect of differential soil drying is reflected more 
strongly on the canopy temperature through differing soil resistances (which is directly related to 
the hydraulic conductivity). Under progressive soil drying a deep rooted crop undergoes a more 
severe stress, although later in cronology. as compared to a shallow rooted crop. The irrigation 
scheduling criteria of 70 percent depletion level of available water means different canopy tempera- 
ture variability for soils with differing texture and crops with differing rooting characteristics. 
While for sandy soils a marked canopy temperature variability may be seen at the 70 percent 
depletion level, very little variability occurs for clayey soils. Alternately, if irrigation scheduling is 
done based on a prescribed canopy temperature variability then the depletion level of the available 
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water would be different for soils with differing texture and crops with differing rooting character- 
istics. Clawson and Btad's ( 1982) criterion for irrigation scheduling (namely, CTV ■ 0.8 K) may not 
be equally effective under all field conditions, particularly in clayey soils. 

In studying the effect of spatial variability of soil hydraulic properties on a bare soil tempera- 
ture, Lascano and van Bavel ( 1 982) found the temperature variabilities to be generally less than 
1 .5 K for a silty clay loam soil. This lead them to conclude that the current technology of infrared 
remote sensing would not be able to detect the effect of soil hydraulic variability. The observations 
of Clawson and Blad (1982) and the present simulation show that canopy temperature variability 
could be in excess of 4K for sandy soils. Such high variability resulting from soil heterogeniety 
should be detectable by the current infrared technology. The present simulation results for the 
clay loam soil would, however, concur with Lascano and van Bavel’s conclusion. A further point 
addressed by Lascano and van Bavel is the accuracy with which the average hydraulic characteristics 
of the soil needs to be known to make model predictions compatible in accuracy with remote 
sensing observations. For their silty clay loam soil, which gave a temperature variability generally 
less than 1 5K. an order of magnitude accuracy in the hydraulic conductivity was suggested. The 
present simulation indicates that the hydraulic conductivity of sandy soils needs to be known with 
a better accuracy than for clayey soils. 

SUMMARY AND CONCLUSIONS 

A soil-plant-atmosphere model for com together with the scaling theory were used to study 
the sensitivity of spatial variation of canopy temperature to soil texture and crop rooting charac- 
teristics. The temperature variability was shown to be more pronounced in a sandy soil as opposed 
to a clayey soil, and in more deeply rooted crop as opposed to a shallow rooted crop. If irrigation 
scheduling is done based on a fixed canopy temperature variability then the level of depletion of 
available water would vary among fields differing in soil texture and crop rooting characteristics. 
Sandy soils would appear to be more amenable to irrigation scheduling via canopy temperature 
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variability as comparedto clayey soils if 70% depletion of available water is recognized as a criteria 
for need of irrigation. 

In clayey soils the spatial variability of canopy temperature would generally be less than 1 K, 
suggesting that infrared radiometric observations with differing spatial resolutions (for example, via 
ground based and air-borne radiometers) should be comparable directly. The spatial variability of 
canopy temperature in sandy soils could be significant enough to require more prudence in com- 
paring radiometric data with differing spatial resolutions. In attempting to 'calibrate* space-borne 
radiometric observations having a coarser resolution by a ground-based observation one needs to 
be cautious about the higher order heterogenieties of soil moisture which could be induced by a 
rainfall event, in addition to the canopy type variations and the presence of bare soils. The effect 
of these higher order (or forced) heterogenieties will be discussed elsewhere. 
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CAPTION TO THE FIGURES 


Figure 1. Schematic illustration of the soil-plant-atmosphere model used in the present simulation. 
The root water uptake is modeled by the van den Honert s (1948) equation; the resist- 
ance analogue model of Federer ( 1 979) is used for changes in plant water storage; 
Monteith’s (1965) equation is used for transpirational flux. The plant water potential 
(0) is calculated as the potential which conserves the water flux in the system. Arrows 
indicate the direction of water flow. Successive steady states are simulated with 
dynamic equations for soil water and plant water storage. 

Figure 2. Comparison of simulated (solid line) and observed diurnal trends of the leaf water poten- 
tial. The observations in (2a) are from Reicosky et al. (1975), and in (2b) from Reicosky 
and Lambert (1978). Concurrently observed microclimatic data are also shown. Plant- 
to-plant differences in the observed potentials are indicated by ‘error’ bars. 

Figure 3. Simulated relationships between leaf and soil water potentials at 1300 hour. The micro- 
climatic data in figure 2b is used in this and all other simulations. The results an for two 
soil textures and two rooting characteristics. Z t ff and L A refer to effective rooting 
depth and root length per unit area, respectively. 

Figure 4. Dependence of the plant water stress index (PWSI. defined in eqn. 10) on the soil water 
potential. The depletion levels of 65 % and 70% of available water are annotated on the 
figures, (a) for clay loam, and (b) for loamy sand. 

Figure 5. The simulated relationship between leaf and soil water potentials at 0600 hour (just 
before sunrise) for the clay loam soil. The meaning of R1 and R2 are clarified in Figure 
3 . A 1 : 1 line is drawn to show that the leaf water potentials are generally lower than the 
soil water potential. 


Figure 6. Simulated leaf stomatal resistances and soil water potentials at 1300 hour for various 
days of the simulation. For the meaning of rooting characteristics, R1 and R2, refer to 
figure 3. 

Figure 7. Simulated relationships of canopy-air temperature differences with soil water potentials 
at 1300 hour. For the meaning of rooting characteristics, R1 and R2, refer to figure 3. 

Figure 8. A simulated range of canopy temperatures at 1300 hour of various days in the simulation. 
A depletion of 70 percent of available water is indicated in the figures. The rooting char- 
acteristics R1 and R2 are clarified in figure 3. 
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